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We investigate the level-density o~(x) and level-spacing distribution p(s) of random matrices M = 
AF 7^ M' where F is a (diagonal) inner-product and A is a random, real symmetric or complex 
Hermitian matrix with independent entries drawn from a probability distribution q(x) with zero 
mean and finite higher moments. Although not Hermitian, the matrix M is self-adjoint with respect 
to F and thus has purely real eigenvalues. We find that the level density (Tf(x) is independent of 
the underlying distribution q(x), is solely characterized by F, and therefore generalizes Wigner's 
semicircle distribution aw(x). We find that the level-spacing distributions p(s) are independent 
of q(x), are dependent upon the inner-product F and whether A is real or complex, and therefore 
generalize the Wigner's surmise for level spacing. Our results suggest f-dependent generalizations 
of the well-known Gaussian Orthogonal Ensemble (GOE) and Gaussian Unitary Ensemble (GUE) 
classes. 



I. INTRODUCTION 

Since their beginning in the field of nuclear physics 
in 1950s, statistical properties of random matrices with 
specific symmetries have been a source of ongoing inves- 
tigations PHI]. Wigner's semicircle law for the eigen- 
value density, and his surmise for the eigenvalue spac- 
ings are some of the most inspired results in random 
matrix theory. For real, symmetric, N x N matrices 
with independent, identically distributed (i.i.d.) entries 
drawn from an arbitrary distribution q(x) that has zero 
mean, variance one, and finite higher moments, Wigner 
showed that the level density in the large-Af limit is 
given by aw(x) = (2/n)y/l — x 2 . He also surmised 
that the level-spacing distribution for random, real, sym- 
metric matrices, p (s) — (ns/2) exp(— 7rs 2 /4), and 
for random, complex, Hermitian matrices, p GUE (s) = 
(32s 2 /-7r 2 ) exp(— 4s 2 /7r), are independent of the underly- 
ing probability distribution q(x). The near-universality 
of these results led to the Gaussian Orthogonal Ensam- 
ble (GOE) conjecture for random, real, symmetric ma- 
trices, the Gaussian Unitary Ensamblc (GUE) conjecture 
for random, complex, Hermitian matrices, and the Gaus- 
sian Symplectic Ensamble (GSE) conjecture for fermionic 
Hamiltonians. The tremendous analytical progress in 
this field is based on these conjectures which imply that 
it is sufficient to consider Gaussian distributed entries 
for the random matrices [4]. We note that the proof of 
Wigner's semicircle law for an arbitrary, non-Gaussian 
distribution q(x) is based on the moment method [T]. 
We also emphasize that the level-spacing distribution 
p (s), although not exact, is an excellent approxima- 
tion to the exact answer obtained for the GOE, and it 
is believed that the result is valid for random, real, sym- 
metric matrices with arbitrary (non-Gaussian) underly- 
ing probability distributions q(x) [4]. 

These ensembles arise from a fundamental constraint: 
the Hermiticity of the random Hamiltonian, that guar- 
antees real eigenvalues, and eigenvectors that are orthog- 
onal with respect to the standard inner-product in quan- 



tum theory [4]- However, a large class of non-Hermitian 
matrices - parity and time-reversal symmetric Hamilto- 
nians [5H7]> rate-equation matrices [8], central potentials 
in momentum space [5], etc. - has real spectra although 
the eigenvectors are not orthogonal under the standard 
inner-product [5]. 

This observation raises the following question: what 
are the properties of random matrices that are self- 
adjoint with respect to a general inner product F = 
and therefore have purely real eigenvalues? A general 
inner-product F, which, by definition, is a positive- 
definite matrix, may represent a system (with non- 
Euclidean geometry) that is not translationally invariant. 
For example, consider a finite disk with rotational sym- 
metry around its center. The disorder potential V(r) at 
a point r on such a disk will only depend on the radial 
distance r from its center, but not on its angular orien- 
tation. Such a disorder can result from sputtering depo- 
sition (for electrons) or from a circular patterned grating 
(for light). The resulting momentum-space Hamiltonian 
for such a potential is not Hermitian; but it has purely 
real eigenvalues !) . Motivated by the generalization of 
this example to D-dimensions, in this paper, we only con- 
sider a diagonal inner product F : Vn x Vn C defined 
by 

JV 

{4>W)F = ^F^ = Y l fi4>i^i (i) 
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where Vn is an A^-dimensional vector space and Fjk = 
djkfj = djkj 0-1 - Note that an arbitrary inner-product 
(•\')f is invariant under transformations T : Vn — > Vn 
that obey FT = F. When F = 1, the standard inner- 
product, this group of transformations corresponds to the 
unitary (orthogonal) group over a complex (real) vector 
space V . Under such an inner-product, an operator M is 
self-adjoint if and only if M = AF where A is Hermitian. 
In this paper, we numerically and analytically investi- 
gate the level density and level-spacing distributions of 
such matrices M, or equivalently M = vF 'AyF = Al\ 



obtained from random, Hermitian matrices A with i.i.d. 
entries drawn from an arbitrary probability distribution 
q(x). 

Our three salient results are as follows: i) The level- 
density <Jp(x) is independent of q{x), is characterized by 
the inner-product F, and therefore, is not the same as 
Wigner's semicircle distribution, ctf(x) 7^ aw(x). ii) For 
bulk eigenvalues, we find that the level-spacing distribu- 
tion is independent of q(x), and is solely characterized by 
F and whether the matrix A is real or complex. Since 
these results are invariant under the group of transfor- 
mations mentioned above, and depend only on whether 
A is real or complex, we use Pp OE (s) and Pp UE (s) to 
represent the level-spacing distributions. Thus Pp OE (s) 
and Pp UE (s) provide -F-dependent generalizations of the 
GOE and GUE universality classes for real and complex 
random matrices respectively, iii) In each case, as D in- 
creases, the level-spacing distribution shifts towards the 
origin and decays slowly compared to the GOE and GUE 
results. 

The plan of the paper is as follows. In Sec. [TT] we 
present numerical data for the level density which show 
a clear F-dependence, and di scus s the properties of the 
level density uf(x). 
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Section 

data for level-spacing distributions in the real and com- 
plex cases. We then present a simple model that quali- 
tatively explains the evolution of level-spacing distribu- 
tions with increasing D > 1 as the inner-product Fp> 
deviates from the standard-inner product, Fd=i = 1- 
We conclude the paper in Sec. |IV| with a brief discus- 
sion. Although the non-Hermitian, random matrices 
M considered in this paper appear similar to those in 
Refs. [IOUl3j . we point out the crucial differences among 
them in Sec. |IV| We emphasize that the probability 
distribution for the random matrices considered here is 
p{M) cx exp(- Tr (A* A)) oc exp(- Tr (F~ 2 M^M))', it is 
not a function of Tr (M^M), and therefore, the tradi- 
tional methods [H QTJ [T2J [T5] used to analyze eigenvalue 
statistics may not be applicable. 



II. LEVEL DENSITY 

We start this section with numerical data that hint at 
our result. For independent entries of the random matrix 
A, we use a Gaussian distribution G(x) with zero mean 
and variance one, or a uniform distribution U(x) with 
zero mean and variance one, or a distribution q$(x) that 
corresponds to random entries Ti = cos OGi + smOUi. 
Note that the distribution qg(x) = q$(—x) interpolates 
continuously from the Gaussian (9 = 0) to the uni- 
form {9 = it/2), has zero mean, variance one, and finite 
higher moments. Since the variance of i.i.d. entries is 
one, the eigenvalue scale for the level density is given by 
A F = 2Tr F/VN (A F = 2V2Tr F/VN) when M is a 
real (complex) matrix. This scale implies that the sec- 
ond moment fj, 2 of the resulting level density is fixed, 

H2 = 1/4 Hang. 
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FIG. 1. (color online) Level-density ap(x) for D — {1, 1.4, 2}. 
The top-left panel reproduces Wigner's result for D = 1, ob- 
tained using N t = 10 3 , N = 3 x 10 3 , N b = 500, and Gaus- 
sian distributed complex entries (0 — 0). The bottom-left 
panel shows the result for D = 1.4 with N t = 500, N = 10 3 , 
Nb — 200, and real entries with 6 — 7r/4. The top-right panel 
shows the result for D = 2, N t = 500, TV = 2 x 10 3 , N b = 200, 
and real entries from a uniform distribution (6 = 7r/2), 
whereas the bottom-right panel shows the result for the same 
D, with N t = 10 3 , N = 3 x 10 3 , N b = 500 and complex 
entries from a distribution q$(x) with 8 = 7r/3. As D in- 
creases, we see that the level-density <tf(x) sharpens near the 
origin, broadens, and, when scaled appropriately, is the same 
for both real and complex random matrices. 



Figure [T] shows the level density for D = {1,1.4,2}. 
We have verified that these results are essentially in- 
dependent of the number of trials Nt 3> 1, the matrix 
size iV > 10 2 , the number of bins 2V&, and the underly- 
ing probability distribution qg{x). The top-left panel in 
Fig. [l] reproduces Wigner's result aw(x) that is expected 
when D = 1. The bottom- left panel (D = 1.4) and the 
top-right panel (D = 2) show that with increasing D, the 
level density up (x) broadens and becomes sharply peaked 
near the origin. The top-right and bottom-right panels 
show that when scaled appropriately, the level density 
ap{x) is the same for real and complex matrices and thus 
depends only on the inner-product F. These numerical 
data strongly suggest that ap(x) is independent of qe(x) 
and is solely characterized by the inner product F. 

To characterize the -D-dependence of the resulting level 
density, we calculate its even moments and compare them 
with the results for the Wigner distribution aw(x) [TlfTp], 
The k th moment (k even) of the level density is given 
by fi k = lim A r^ 00 ( Tr ({AF) k )/N) A where (• • • ) A denotes 
random averaging. It is straightforward, but tedious, to 
calculate the non-vanishing averages systematically; for 
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FIG. 2. (color online) Even moments fik of the level density 
ctf(x) scaled by their values fikw for the Wigner distribution, 
for D = {1, 1.5, 2}. Each numerical data-point (solid symbols) 
is obtained by using a different set of values for (Nt, N, qg(x)). 
As D increases the ratio Hk/^kw ~ I)( fe / 2_1 ) increases, con- 
sistent with the broadening of the level density (Jf(x) relative 
to the D = 1 result. 
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Note that when F = 1 the sum of terms in a square 
bracket gives the corresponding Catalan number [TJ 116) . 
Since F = diagO' 23 " 1 ), TrF m - Arm(£>-i)+l for D > 
and finite moments require D > 1. Fig. [2] shows the nu- 
merically obtained even moments scaled by their values 
for D = 1, Hu/p^kW (solid symbols), and the analytical 
results (open symbols) for D = {1,1.5,2}. The mono- 
tonic increase in the even moments with D is consistent 
with the broadening of the level density ap(x) relative to 
the Wigner 's result aw(x). Although an exact solution 
for ap(x) is known through implicit equations |17j . our 
(numerical and analytical) approach provides detailed in- 
formation. 

We recall that the level-density uf(x) is independent 
of the underlying probablity distribution q(x) only when 
q(x) has finite higher moments; for example, even for 
F = 1, if q(x) is the Cauchy distribution, the result- 
ing level density is not the same as the Wigner distri- 
bution <7w{x) [16,. We also point out that the exis- 
tence of the asymptotic limit in Eqs.pJ-Q is determined 
by the large-A^ behavior of the inner-product elements 



Fjk = Sjkfj- In particular, if liniAr^oo TrF is finite or 
diverges as N a with a < 1, it follows that the level den- 
sity ap(x) has divergent higher moments and therefore it 
cannot have a compact support. Numerical data obained 
by using F jk = Sjuj ' 1 with D < 1 or F jk = 5jke~ oa 
show that the support of the level density <Jf,n{x) widens 
as N increases. Thus, universal results for the level den- 
sity <jf{x) are obtained only when the underlying proba- 
bility distribution q{x) and the inner-product F obey the 
constraints established here. 



III. LEVEL-SPACING DISTRIBUTION 



by 



The level spacing distribution pf(s) is formally given 



Pf(s) = lim 




(•5) 



where Xi(M — AF) are the eigenvalues listed in as- 
cending order, > A,, and AA avg = Ylj=i OVkl - 
Xj)/(N — 1) is the average spacing between the eigenval- 
ues. Note that the eigenvalue scale Ap does not explicitly 
enter Eq.j5J). 

Figure [STshows numerically obtained level-spacing dis- 
tribution for the bulk (> 10%) eigenvalues of random, 
real, N x N matrices M = AF for D = {1, 2.5, 3}. Note 
that the level-spacing scale is chosen so that the average 
level-spacing is one: J °° spp(s)ds = 1. We have verified 
that these results are independent of Nt 3> 1, N > 10 2 , 
and qg(x) used to create the matrix A. When D = 1, 
F is the identity matrix and the numerical results (black 
squares) match well with Wigner surmise p GOE (s) (black 
solid line). As D increases, we see that the level-spacing 
distribution becomes broader, the weight of the distribu- 
tion near the origin increases and so does its slope at the 
origin. We recall that this distribution is invariant under 
real matrix transformations O that obey O t FO = F. 
Since when F = 1, they correspond to the group of or- 
thogonal matrices, we use the notation Pp OE (s) to de- 
note the level spacing results that are invariant under 
the F-dependent group of real matrix transformations. 

Figure [4] shows the level-spacing distribution for the 
bulk (> 10%) eigenvalues of random, complex, N x N 
matrices M = AF for D — {1,2.5,3.0}. When D = 1, 
the numerical results (black squares) match with the 
Wigner's surmise p (s) (black solid line). As D in- 
creases, we see that the level-spacing distribution be- 
comes broader, the weight of the distribution near the ori- 
gin increases and so does its second derivative at the ori- 
gin. This distribution is invariant under complex matrix 
transformations U that obey WFU — F. Since when 
F = 1, they correspond to the group of unitary matrices, 
we use the notation Pp UE (s) to denote the level spacing 
results that are invariant under the _F-dependent group 
of complex matrix transformations. These numerical re- 
sults strongly suggest that the level-spacing distributions 
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FIG. 3. (color online) Level-spacing distribution Pp OE (s) for 
random, real, N x N matrices M = AF with D = {1, 2.5, 3}. 
The D = 1 results (black open squares) are obtained by using 
N t = 10 3 , N = 2 x 10 3 , N b = 500, and a Gaussian distri- 
bution, and match well with Wigner surmise p GOE (s) (black 
solid line). The D — 2.5 results (red circles), obtained by 
using N t = 10 3 , N = 3 x 10 3 , N b = 500, and 9 = ir/4, show 
that the level-spacing distribution maximum shifts towards 
the origin, its slope increases, and it broadens. The D — 3 
results (blue stars) are obtained by using Nt = 100, N = 10 4 , 
Nb = 500, and a uniform distribution. 
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FIG. 4. (color online) Level-spacing distribution pp UE (s) 
for random, complex, N x N matrices M — AF with D = 
{1,2.5,3}. The D = 1 results (black open squares) are ob- 
tained by using Nt = 10 3 = N, Nb = 200, and a Gaus- 
sian distribution, and match well with the Wigner's surmise 
p GUE {s) (black solid line). The D = 2.5 results (red circles), 
obtained using JV* = 10 3 = N, N b = 500, and 9 = 7r/3, show 
that the level-spacing distribution maximum shifts towards 
the origin, its second derivative at the origin increases, and 
it broadens. The D — 3 results (blue stars) are obtained by 
using Nt = 100, N = 3 x 10 3 , N b = 500, and a uniform 
distribution. 



p GOE (s) and p GUE {s) are independent of the underlying 
probability distribution qg(x) and are F-dependent gen- 
eralizations of the GOE and GUE universality classes. 

The D-dependence of the large- N level-spacing distri- 
butions may be qualitatively understood by exploring it 
for 2 x 2 matrices M = AF. Here, A is a 2 x 2 ran- 
dom, real symmetric or complex Hermitian matrix with 
entries drawn from a Gaussian distribution with zero 
mean, variance a, and F = diag(l, fa) where we choose 
/2 = N ^ 1 ^> 1 to mimic the effects of a non-standard 
inner-product F. We remind the Reader that when 
D = 1, this procedure gives excellent approximations, 
p GOE {s) and p GUE (s), to the large-iV results [3 HH]. 
When D > 1, for real, 2x2 matrices M = AF, the 
result is [TH] 

P^°2 E (s) = - RS e^ 2 I Q {B^) (6) 
ap 

where a 2 = 16a 2 f 2 , /3 2 = 8<r 2 (l + f 2 ), B ± = (a- 2 ±(3~ 2 ), 
and Io(x) is the modified Bessel function of first kind. 
The variance a is chosen so that the mean level-spacing 
is unity. In the limit f 2 = N D - X » 1, B± -> vr/ 2 /32, 
(B + — _B_) —> 1/ir, and the constraint on the mean-level 
spacing leads to p G ° E { s ) — ( n \/ f2/32)sexp(—s 2 /n). 
Thus, as D increases, the distribution p GGE {s) shifts to- 



wards the origin, the slope of the distribution near the 
origin diverges, and at large s, the distribution decays 
more slowly than p GOE ~ exp(— 7rs 2 /4). When the 2x2 
matrix A is complex, the result for the level-spacing dis- 
tribution is given by 

P2X2 (s) = \ -/X i/s e 7 7) 

V 7T ST 

where = 4cr 2 / 2 , i' -2 = 2ct 2 (1 + / 2 2 ), D(x) is the 
Dawson's integral which satisfies D(x) ~ x when i«l 
and D(a;) — l/(2x) when x > 1, and r 2 = (/i 2 - ^ 2 )/2. 
When /2 = N ^ 1 3> 1, the c onst raint on the mean- 
level spacing implies that cr/2 \/4/7r = 1. This leads to 

P2x2 E {s) = v / 8/2/7r 3 se- s2/7r L>(s v / / 2 /27r). Therefore, as 
D > 1 increases, the distribution p G ^ E {s) shifts towards 
the origin, the second derivative of the distribution at 
the origin diverges, and at large s the distribution decays 
more slowly than p GUE ~ exp(— 4s 2 /-7r). 

Figure [5] shows the level-spacing distributions p G ® E {s) 
(blue squares thick line), the corresponding Wigner sur- 
mise p GOE (s) — (tts/2) exp(— 7rs 2 /4) (blue solid thin 
line), the distribution for complex matrices p G ^ E { s ) 
(red dotted line), and the corresponding Wigner surmise 
p GUE (s) = (32s 2 /7r 2 )exp(-4s 2 /7r) (red dashed line) ob- 
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FIG. 5. (color online) Normalized level-spacing distribution 
P?x2 B > Eq.ljfj]), (blue squares thick line), the Wigner surmise 
p GOE ^ blue golid thin line ^ the distribution p^xf , Eq.jT]), 
(red dotted line), and the Wigner surmise p GUE (red dashed 
line). In each case, as fa 3> 1 increases, the distribution P2x2 B 
(P2x2 ) shifts towards the origin, its slope (second derivative) 
at the origin increases, and the distribution broadens. 



tained by using /2 = 100 ^> 1. We emphasize that al- 
though the 2 x 2-case analytical results qualitatively corre- 
spond to the large- N numerical data in regions s <C 1 and 
s > 1, they are not a good approximation in the inter- 
mediate range, and the analytical answer for the large- N 
level-spacing distributions Pp OE (s) and Pp UE (s) remains 
unknown. 



IV. DISCUSSION 

In this paper, we have investigated the level-density 
and level-spacing distributions for matrices M — AF ^ 
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We have shown that the level density ctf(x) is depen- 
dent upon, and solely characterized by the inner-product 
F. With increasing D, the level density o-f(x) develops a 
peak at the origin, and broadens its base, but continues to 
have a compact support. We found that the level-spacing 
distribution Pf(s) for such matrices is dependent upon 
the inner-product F, and upon whether the matrix A 
is real or complex. As D increases, the level-spacing dis- 
tributions p ( f OE {s) and p'f UE {s) shift towards the origin, 
become steeper near the origin, and broaden. Our numer- 
ical results and qualitative analysis thus strongly suggest 
F-dependent generalizations of the GOE and GUE uni- 
versality classes; in each case, however, analytical expres- 
sions for the independent n-point correlation functions 
remain unknown [4]. 

The analysis in this paper is limited to a specific func- 
tional form of the inner-product F. Our results, how- 
ever, raise the broader question of universality classes 
for random matrices that are self-adjoint with respect 
to a generic, non-diagonal, positive-definite, inner prod- 
uct F. The random, non-Hermitian matrices M = AF 
are similar to those in Refs. [101 [11] but for the crucial 
difference that, in contrast to Refs. (TUJITT], the eigen- 
values of M = AF are always purely real. The probabil- 
ity distribution p(M) oc exp(— Tr (A* A)) is also similar 
to the measure for random, generalized Wishart matri- 
ces [THEI]. I n contrast to Refs. [HI [T3], however, we 
have explored the eigenvalues of the matrix M = AF, 
and not the positive-definite matrix A^A. The gener- 
alizations of Wigner's semicircle law and the GOE and 
GUE surmises for level-spacing distributions, presented 
in this paper, suggest that the statistical properties of 
random, self-adjoint matrices, although independent of 
the underlying probability distribution q(x), are richly 
varied. 



ACKNOWLEDGMENTS 

W.A.K. was supported by a IUPUI Undergraduate Re- 
search Opportunities Program (UROP) grant. 



[8] C. Timm, Phys. Rev. E 80, 021140 (2009). 

[9] See Eqs.(10) and (11) in W.A. Karr, C.R. Jamell, and 

Y.N. Joglekar, Am. J. Phys. 78, 407 (2010). 
[10] P. Shukla, Phys. Rev. Lett. 87, 194102 (2001). 
[11] A.M. Garcfa-Garci'a, S.M. Nishigaki, and J.J.M. Ver- 

baarschot, Phys. Rev. E 66, 016132 (2002). 
[12] A.Y. Abul-Magd, G. Akermann, and P. Vivo, J. Phys. A 

42, 175207 (2009). 
[13] Vinayak and A. Pandey, Phys. Rev. E 81, 036202 (2010). 
[14] L.A. Pastur, Lett, in Math. Phys. 25, 259 (1992). 
[15] E Brezin, V. Kazakov, D. Serban, P. Wiegmann, and A. 



6 



Zabrodin, Applications of Random Matrices tn Physics 
(Springer, New York, 2006). 
[16] S.J. Miller and R. Takloo-Bighash, An Invitation to Mod- 
ern Number Theory (Princeton University Press, Prince- 
ton, 2006). 

[17] Z.D. Bai and L.X. Zhang, J. Multivariate Anal. 101, 1927 



(2010). 

[18] D. Herman et al, Phys. Rev. B 76, 195448 (2007). 
[19] P. C. Huu-Tai, N.A. Smirnova, and P. Van Isaker, J. 
Phys. A 35, L199 (2002). 



